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We present a systematic study of twist-3 light-cone distribution amplitudes of vector 
mesons in QCD, which is based on conformal expansion. A complete set of distribution 
amplitudes is constructed for p, u, K* and <fi mesons, which satisfies all (exact) equations 
of motion and constraints from conformal expansion. 



Amplitudes of hard (light-cone dominated) exclusive processes in QCD are expressed 
by the factorization formulae, which separate the short- distance dynamics from the long- 
distance one|I|]. For those processes producing a (light) vector meson V (= p, u, K*, 0) in 
the final state, like exclusive semileptonic or radiative B decays (B — > Vev, B — > V + 7) 
and the hard electroproduction of vector mesons (7* + N — > V + N'), a long-distance part 
involving the final vector meson is given by the vacuum-to-meson matrix element of the 
nonlocal light-cone operators (0\if}(z)r\ l 'if)(—z)\V), corresponding to universal nonpertur- 
bative quantity called the light-cone distribution amplitudes (DAs). {z p is a light-like 
vector z 2 = 0, and A 4 and T denote various flavor and Dirac matrices. We do not show 
the gauge phase factor connecting the quark and the antiquark fields.) 
/\! The DAs of higher twist are essential for systematic study of preasymptotic corrections 

to hard exclusive amplitudes, and are also interesting theoretically because they contain 
new and direct information on hadron structure and the dynamics of QCD. 

Here we will explicitly consider the "chiral-odd" DAs of the charged p meson with 
momentum P p and polarization vector e^' (P 2 = mi, e^ ■ e^ x ' = —1, P ■ e^) = 0). 
The "chiral-even" DAs can be treated similarly 0; and extension to other vector mesons 
is straightforward (see below). The relevant quark-antiquark DAs are defined with the 
chirality-violating Dirac matrix structures V = {a^, 1} as 
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\u(z)d(-z)\p-(P,X)) = -1 fj - f P ^^)ml(e^ . z ) [ du e ^h\°\u,p 2 ). (2) 



In ([!]), we neglect the Lorentz structures corresponding to the twist-4 terms for simplic- 
ity (see Ref.0 for the complete expressions). Our Lorentz frame is chosen as P ■ z = 
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P + z . The nonlocal operators on the l.h.s. is renormalized at scale fi. We set £ = 
u — (1 — u) = 2u — 1, and / p and /J are the usual vector and tensor decay constants as 
(0|«(0) 7 ^(0)|p-(P, A)) = f p m p eW and (0|u(0)<r^d(0)|p-(P, A)) = ifJ(e^P u - e^P,). 

We wrote e^ = ejjj + eg* with ejj = [P^ - z p m 2 p /(P ■ z)](e< A > ■ z)/{P ■ z), so that the 
DAs with the subscripts || and _L describe longitudinally and transversely polarized p 
mesons, respectively. (j)± is of twist-2, while hi and hf are of twist-3. These DAs are 
dimensionless functions and describe the probability amplitudes to find the p in a state 
with quark and antiquark, which carry the light-cone momentum fractions u and 1 — u, 
respectively, and have a small transverse separation of order 1/p. 

The quark-antiquark-gluon DAs can be defined similarly; they are of twist-3 and higher. 
There exists one chiral-odd DA of twist-3, which is given by 

(0\u{z)a^z v gG m (vz)z^d{-z)\p-(P,X)) 

= {P ■ z)(e (A) ■ z)fl p m p f Va e- iPz(a "- Qd+1,Q9) T(a,/i 2 ), (3) 

where G m is the gluon field strength tensor, Va = da.dda u da g 5{l — otd — ot u — a g ), and a 
denotes the set of three light-cone momentum fractions: a^ (quark), a u (antiquark), and 
a g (gluon). f£ is the three-body tensor decay constant 0, so that T is dimensionless and 
is conveniently normalized as JT>a(ad — a u )T(a) = 1. The three-particle DA T describes 
a higher component in the Fock-state wave function with additional gluon. 

The basis of DAs defined above is overcomplete due to the constraints from the QCD 
equations of motion. Using [%Jp — m q )q(x) = (g = u, d), we obtain 

Id r 1 

-a;"- — u(x) [nfn,nf v ] ± d(-x) = / dvv ( ^ lTl)/2 u(x)a flu x u gG m (vx)x ri d(-x) 

x v d 



2 % 



[u{x + y) [fn, 7j T d(-x + y)} _^ + i(m u ± m d )u(x)^d(-x),(4) 



where [7^, ^ v ] ± = 7 M 7 I/ ±7 !/ 7 At . In the light-cone limit x 2 -^ 0, the vacuum-to-meson matrix 
elements of (||) yield a system of integral equations between two- and three-particle DAs 
defined above. We note that the total derivative term induces mixing between hi and 
hi? . We can solve0 these coupled integral equations in a form (i = s,t) 

hf{u,p 2 ) = Cdv K$ w (u,v)<f>_ L (v,n 2 ) + Cva K®(u,a)T(a, /1 2 ), (5) 
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where K^ (X = WW, g) are independent of /1. We omit the quark mass correction 
term for simplicity. The first term on the r.h.s. is the twist-2 contribution and thus the 
analogue of the Wandzura-Wilczek piece of the nucleon structure function g2(x, Q 2 ). 

To analyze ([5]), we introduce the conformal partial wave expansion in the m q — > 
limit [0,^J . The conformal expansion of light-cone DAs is analogous to the partial wave 
expansion of wave functions in standard quantum mechanics (QM). In conformal expan- 
sion, the invariance of massless QCD under conformal transformations substitutes the 
rotational symmetry in QM. In QM, the purpose of partial wave decomposition is to sep- 
arate angular degrees of freedom from radial ones (for spherically symmetric potentials). 



All dependence on the angular coordinates is included in spherical harmonics which form 
an irreducible representation of the group 0(3), and the dependence on the single remain- 
ing radial coordinate is governed by a one- dimensional Schrodinger equation. Similarly, 
the conformal expansion of DAs in QCD aims to separate longitudinal degrees of freedom 
from transverse ones. All dependence on the longitudinal momentum fractions is included 
in terms of certain orthogonal polynomials which form irreducible representations of the 
so-called collinear subgroup of the conformal group, SL(2,R). The transverse-momentum 
dependence (the scale-dependence) is governed by simple renormalization group equations: 
the different partial waves, labeled by different "conformal spins" , behave independently 
and do not mix with each other. Since the conformal invariance of QCD is broken by 
quantum corrections, mixing of different terms of the conformal expansion is only absent 
to leading logarithmic accuracy. Still, conformal spin is a good quantum number in hard 
processes, up to small corrections of order a 2 . 

The conformal expansion of the DAs on the r.h.s. of (0) reads 
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where u = 1 — u, £ = 2m — 1, a g = l — ad — a u , and C^ 2 and J k i are particular Gegenbauer 
and Appell polynomials [0. Using orthogonality relations of these polynomials, the expan- 
sion coefficients are expressed by matrix elements of local conformal operators: a n and uj k i 
are given by local operators of conformal spin j = n + 2 and j = k + I + 7/2, respectively 
(three-particle conformal representations are degenerate). Thanks to conformal symme- 
try, Kx °f © are resolved into the superposition of the terms of definite conformal spin. 
As a result, h» are given by the conformal expansion, where each expansion coefficient 
is expressed in terms of a n and u k i with the corresponding spin j. 

Now all DAs are expressed, order by order in the conformal expansion, by indepen- 
dent matrix elements a n and u) k i- From (|5p and (^), the /^-dependence of the DAs is 
governed by that of a n (« 2 ) and ^/(m 2 ). This is determined by the renormalization of 
the corresponding local conformal operators, and is worked out in the leading logarithmic 
approximation PQ. The results indicate that the relevant anomalous dimensions increase 
as ~ In j. This means that only the first few conformal partial waves contribute at suf- 
ficiently large scales. Therefore, the truncation of the conformal expansion at some low 
order provides a useful and consistent approximation of the full DAs. 

We take into account the partial waves with j < 9/2, where the terms with n < 2 and 
k + 1 < 1 are retained in the first and second equations of (0). Correspondingly, we getp| 
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+ 36 + (3uu + u\nu + uliau) + 3<5_(tiln-u — uliau), (7) 

hf(u) = 3e 2 + ^(3e 2 -l) + ^ 2 (5e 2 -3) + ^C 3 (3-30e 2 + 35e 4 ) 
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+ -5 + (l + £\nu/u) + -5^(2 + \nu + \nu). (8) 

Here ^3 = fl p /{.fj rn p)^ an d we incorporated the quark mass corrections proportional 
to 5± = f p {m u ± md)/{fjmp). The results for other vector mesons u, K* and 4> are 





Figure 1. Two-particle twist-3 DAs for the p, K* and <\> mesons: (a) hi (w) of (|7|); (b) 

hi [u) of (§). "p (WW)" denotes the 'Wandzura-Wilczek" type contribution given by 
the first three terms of (^) and (§) for the case of the p meson. 



obtained by trivial substitutions. We emphasize that these results provide a consistent 
set of DAs involving a minimum number of independent nonperturbative parameters. 
These parameters are calculated from QCD sum rules taking into account SU(3)-breaking 
effects[Q. The resulting DAs are plotted in Fig. |l| The comparison between the solid and 
dotted curves shows that the three-particle contributions (the term with ( 3 in (0), (|[)) 
are important and broaden the distributions. On the other hand, the quark mass effects 
are not so large except for the end-points u — > and u — > 1 for h\ (u) . 

In conclusion, we have developed a powerful framework for hard exclusive processes, 
which allows to express the higher twist DAs in a model-independent way by a minimum 
number of nonperturbative parameters. Combined with estimates of the nonperturbative 
parameters, this leads to model building of the DAs consistent with all QCD constraints. 
Our formalism is applicable to arbitrary twist and other light mesons||, and our results 
are immediately applicable to a range of phenomenologically interesting processes |7[]. 
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